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ABSTRACT 

We analyze the variability in accretion-related emission lines for 40 Classical T Tauri stars to probe 
the extent of accretion variations in young stellar objects. Our analysis is based on multi-epoch high- 
resolution spectra for young stars in Taurus- Auriga and Chamaeleon I. For all stars, we obtain typically 
four spectra, covering timescales from hours to months. As proxies for the accretion rate, we use the 
Ha 10% width and the Ca II-A8662 line flux. We find that while the two quantities are correlated, 
their variability amplitude is not. Converted to accretion rates, the Ca II fluxes indicate typical 
accretion rate changes of 0.35 dex, with 32% exceeding 0.5 dex, while Ha 10% width suggests changes 
of 0.65 dex, with 66% exceeding 0.5 dex. We conclude that Ca II fluxes are a more robust quantitative 
indicator of accretion than Ha 10% width, and that intrinsic accretion rate changes typically do not 
exceed 0.5 dex on timescales of days to months. The maximum extent of the variability is reached 
after a few days, suggesting that rotation is the dominant cause of variability. We see a decline of the 
inferred accretion rates towards later spectral types, reflecting the M vs. M relationship. There is 
a gap between accretors and non-accretors, pointing to a rapid shutdown of accretion. We conclude 
that the ~ 2 orders of magnitude scatter in the M vs. M relationship is dominated by object-to-object 
scatter instead of intrinsic source variability. 

Subject headings: accretion, accretion disks — circumstellar matter — stars: formation — stars: 
low-mass, brown dwarfs — planetary systems 



1. INTRODUCTION 

Young low-mass stars that are still accreting while con- 
tracting towards the main-sequence are identified ob- 
servationally with T Tauri stars (e.g. iBertoutI Il989l : 
lAppenzeller fc: Mundtl Il989f ). One of th e ir defining 
properties is variab ility (e.g. iKnottI 118911 : iJovl Il945t 
iRvdgren et al.l fl97^ . which can be traced to changes 
in continuum excess flux caused by varying accretio n 
flow from a circumstellar disk (|Hartigan et al.l |1991D : 
other possible sources of variability include cool spot 
rotation and extinction events. In addition to the 
continuum, this variability is seen in a number of 
emission lines, including the Hydrogen Balmer series, 
and has been attri buted to rotation, variable accre- 
tion, a nd winds (e.g . Johns & Basri 1995: 'Batalha et alj 
120011: lAlencar fc Ba talha 2002; Alcncar ct al. 200^ 
Photometric monitoring campaigns for large num- 
bers of objects have revealed that hot spots formed 
by gas accretion, co-rotating with the objects, are 
one of the mos t important sources of var iability in 
young stars fe.g. iHerbst et al.lll994HBouviere t al. 199;!; 
iFernandez fc Eiroalll99ir However, it remains unclear 
to what extent the photometric and spectroscopic vari- 
ability seen in typical T Tauri stars is directly related to 
changes in the accretion rate. As we argue below, this 
issue can be addressed by spectroscopic monitoring, but 
so far spectroscopic studies have been limited to small 
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samples. 

Inferences about accretion are usually made in the 
framework of the magnetospheric accretion scenario (see 
the review bv lBouvier et al.l[2007f ). Recent observations 
have shown the presence of a correlation between in- 
ferred mass accretion rate and central object mass, ex- 
tending over se veral orders of magn i tude: M oc M" , 
with a ~ 2 (e.g iMuzerolle et al.l[2003l : iNatta et al.ll2004l : 
iMohantv et al.|[2005f) . Apart from this correlation, all M 
vs. M plots feature a large scatter: at any given object 
mass, the accretion rates show a dispersion of about two 
orders of magnitude. 

A number of ideas have been put forward to 
interpret these findings. The accretion rate vs. 
mass correlation has been attributed to a Bondi- 
Hoyle flow to the star-disk system (iPadoan et al.l 
|2005|). to initial rotatio nal velocities of collapsing cores 
(iDullemond et al.]l2006D. to a dispe rsion in disk param- 
eters () Alexander fc Armitcigil2006D, to a complex mag- 
netic field geometrv ("Grego rv et al.l 12006V an d to a de- 
clining disk ionization with stellar mass (Muz crolle et al.] 
l2003HMohantv et al.ll2005^ . Alternatively, it has been 
argued that the correlation is not physical, but rather 
reflects selection effe cts; therefore, the domi nant feature 
is the scatter itself (jClarke fc Pringld[200l l. AU these 
scenarios predict specific values for a (0-2) and specific 
properties of the scatter. 

A basic question about the M vs. M relation is whether 
scatter around it is due to object-to-object variations 
(i.e. the accretion rate for any individual object re- 
mains mostly constant) or due to a variable acc retion 
rate in any given object (jScholz fc Javawardhanal [2006'). 
At least for some objects, strong changes in accretion 
related lines have been reported, e.g. for the brown 
dwarf 2MASSW J1207334-393254 (|Scholz et all [20051 : 
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IStelzer et al.ll2007l ). but it is not clear if such variations 
are common or not. Since most previous studies were 
based on single-epoch measurements of accretion rates, 
it was impossible thus far to address this issue. 

Previous work has established that a number of op- 
tical emission lines originate in the accretion flow and 
are affected by the accretion rate. In particular, the 
Ha 10% width and the Ca II-A8662 lin e flux are found 
to correlate well with the accretion rate ([MuzeroUe et al.l 
ll998HNatta et al.ll2004l : rMohantv et al.ll2005f) . These em- 
pirical indicators facilitate studies of accretion, and for 
the first time, allow us to investigate accretion in multi- 
epoch spectra for large samples of objects. Here, we use 
both diagnostics to probe the intrinsic accretion variabil- 
ity in young stars in Taurus- Auriga and Chamaeleon I, to 
provide new observational limits for the aforementioned 
scenarios. 

2. OBSERVATIONS AND DATA ANALYSIS 

We obtained multi-epoch high-resolution optical spec- 
tra of 40 members in the ^ 2 Myr old Chamaeleon I and 
Taurus-Auriga star forming regions. The targets span 
the spectral types from F2 to M5 based on published 
classifications, and consist of accretors wit hout suspected 
close companions selected from the list of 'Nguve n et al.l 
(j2009). The data were collected using the echelle spec- 
trograph MIKE (jBernstein et al.ll2003l ) on the Magellan 
Clay 6.5 meter telescope at the Las Campanas Observa- 
tory, Chile on 15 nights during four observing runs be- 
tween 2006 February and 2006 December. Each target 
was observed typically at four epochs with baselines of 
hours, days, and months (only one target was observed 
just twice; 24, 12 and 3 targets were observed four, five 
and six times, respectively). For the data reduction, we 
used customized routines running in the ESO-MIDAS 
environment (described in detail in Brandeker et al., in 
preparation). 

The Ha 10% widths were determined as follows. First, 
we estimated the continuum level by linearly interpo- 
lating flux values in the range of 500 km s~^ to 1000 
km s~^ on either side of the line. Next, the maximum 
flux level of Ha emission was measured with respect to 
this continuum level. Finally, the crossing points of the 
Ha emission with the 10% flux level were identified, and 
the width was measured. Sometimes absorption compo- 
nents in the Ha emission line profile falls below the 10% 
flux level. In such cases, we ensured that we measured 
the width consistently for all epochs, i.e. we measured 
widths to the edges of any blue- or red-shifted absorption 
features. We did not correct for underlying photospheric 
absorption in Ha. 

We derived the Ca II-A8662 emission fluxes (.Fcan) 
from the observed emission equivalent widths. To deter- 
mine the widths, we integrated the emission above the 
continuum level. For emission profiles attenuated by a 
broad absorption feature, we used the median flux within 
0.2 A of the absorption minima as an approximate con- 
tin uum level for in t egrati on, similar to what was done 
bv iMuzerolle et all (|1998l ). To infer the emission fluxes 
from the equivalent widths, we must know the underly- 
ing photospheric continuum flux. We used the contin- 
uum flux predicted by the PHOENIX synthetic spectra 
for a specified Tog and surface gravity. We inferred Tog 
from our spectral types, and assumed a surface gravity 



of log g = 4.0 (cgs units). In our estimate, we ignore veil- 
ing, which could lead to an under estimate of the line flux . 
Indeed, for five targets shared bv lMohantv et al.l (|2005D . 
our results were lower by 0.05 to 0.41 dex. We tried to 
measure veiling from our spectra, but we found the rel- 
atively poor S/N prevented us from reaching sufhcient 
accuracy (S/N of our spectra was typica lly ~ 20 at Ha, 
where as literature studies for veiling, e.g. lHartigan et all 
()1989t ). use spectra with S/N > 75). However, for our 
purposes of studying variability, the bias due to veiling 
is not important. 

3. ACCRETION INDICATORS AND THEIR VARIABILITY 

Previous studies have shown that the Ha 10% width 
and the Ca II-A8662 flux are correlated with accretion 
rates determine d using the tradit i onal method based o n 
optical veiling (jNatta et aLll2004HMohantv et al.ll2005D . 
wit h Ca II showing significantly less scatter around the 
fit (|Herczeg fc Hillenbrandll200l . In Fig. [J we compare 
these two accretion indicators for our sample, with the 
points set to the average values from the multi-epoch 
spectra, and the 'error bars' indicating the range of val- 
ues (the ranges are dominated by variability, not by mea- 
surement uncertainty). We find a clear linear correlation 
between the Ha 10% width and log flux in Ca II-A8662, 
with the deviations comparable to the scatter. This pro- 
vides reassurance that both parameters are mainly deter- 
mined by the same physical quantity, which we identify 
as the mass accretion rate based on published findings. 
For individual objects, the degree of correlation between 
the two accretion indicators varies, which may be ex- 
plained by the varying extent to which accretion affects 
different emission lines. We do not see a significant corre- 
lation between the variability in the accretion indicators 
and stellar mass, Spitzer IR excess, or star forming region 
in our sample. 

In Fig. [21 we compare the variability in the two ac- 
cretion indicators. Here we use (max — min) for the 
Ha 10% width and (max/min) for the Ca II-A8662 flux, 
because these quantit ies are directly proportional to ac- 
cretio n rate changes ()Natta et al.l 1200^ : iMohantv et al.l 
[2005h . One sees that the data are not correlated (cor- 
relation coefficient ~ 0.02), indicating that at least 
one of these parameters does not reflect changes in the 
accretion rate. 

Using the correlations logAf = 1.061ogJ^can — 15.40 
of IMohantv et all ([20051), and logM = -12.89 + 
0.0097Hal0% [km s'^] of lNatta et aH (|2004D . we con- 
vert the variability seen in the spectroscopic indica- 
tors to accretion rate changes. Based on Ca II-A8662 
fluxes, the median accretion rate variability in our spec- 
tra is 0.35 dex, with 13/40 or 32% exceeding 0.5 dex and 
only 1/40 or 2.5% exceeding one order of magnitude. 
These numbers should be treated as upper limits, as 
the uncertainties in measuring Ca II fluxes likely con- 
tribute ^ 0.1 dex to the scatter. In contrast, using the 
Ha 10% width gives a median variability of 0.65 dex, with 
the majority exceeding 0.5 dex (26/40 or 65%), and still 
a substantial fraction exceeding one order of magnitude 
(16/40 or 40%). 

This result indicates clearly that the Ca II-A8662 flux 
is a more robust quantitative diagnostic of accretion rate. 
Our empirical result is c onsistent with t he m odels of 
Classical T Tauri stars of lAzevedo et al.l (|2006l ). which 
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Fig. 1. — Ca II-A8662 fluxes vs. Ha 10% widths in our sam- 
ple: shown are the averages from our multi-epoch spectra with the 
error bars indicating the minimum-maximum range. From linear 
regression analysis, the probability that the data are not correlated 
is < 10~®; the best linear fit is overdrawn. 
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Fig. 2. — Variability in spectroscopic diagnostics of accre- 
tion measured in our multi-epoch spectra for individual objects: 
(max/min) in Ca II-A8662 fiux vs. (max — min) of the 10% width. 
Both quantities are supposed to be directly related to accretion rate 
changes. We do not find a significant correlation in this dataset 
(correlation coeflicient ~ 0.02). 



find that Ca II broad components are formed predom- 
inantly in the accretion flow, and track accretion rate, 
while the hydrogen lines are much more affected by stel- 
lar winds. 

The main reason to use Ha for accretion rate mea- 
surements is that it is much easier to measure, and is 
sensitive to lower levels of accretion. In comparison with 



line fluxes or equivalent widths, the 10% width also has 
the advantage of not being heavily affected by uncer- 
tainties in estimating the underlying continuum. This 
property is beneficial especially for faint objects: the typ- 
ical variation in 10% width from continuum uncertainty 
is ~ 15 km s~^. Based on our findings, however, the 
10% width should not be trusted for quantitative mea- 
surements. Nevertheless, it is still useful as a qualitative 
indicator of accretion. 

A major problem in the 10% width measurement is 
its strong dependence on the line profile. The major- 
ity of the accreting stars in our sample show absorption 
components in their Ha profile, caused by relatively cool 
gas, e.g. in stellar or disk winds, or parts of the accretion 
flow, seen in projection against the hot shock front. Since 
these absorption features are variable in intensity and in 
position relative to the emission, the maximum intensity 
and therefore the 10% level can vary considerably from 
epoch to epoch. This effect can cause large changes in 
10% width not necessarily associated with accretion rate 
changes. 

To probe the timescales of variability in the two line 
indicators, we used our time series baselines of hours 
to months: we examined indicator changes between all 
combinations of epochs for each object, e.g. if an ob- 
ject was observed at epochs A, B & C, then we re- 
viewed changes between epochs A-B, A-C, and B-C. We 
show the changes graphically in Fig. [31 Both for the 
Ha 10% width and for the Ca I1-A8662 fiux, we obtain 
consistent results: the extent of the variability increases 
on timescales ranging from hours to several days. On 
longer timescales, the amount of variability saturates. 
Thus, the dominant timescale for accretion- related vari- 
ability is several days. This idea suggests the variabil- 
ity is determined near the star, where the timescales are 
sufficiently short. Probably, one major factor is rota- 
tion, since typical rotati on periods in youn g stars are 
in the range of 1-10 d (|Herbst et all |2007[ ). In addi- 
tion, there may also be a connection to the character- 
istic infall timescale, typically hours, or magnetic field 
reconnection events, perhaps several days. Note that 
we cannot probe long-term variations on timescales of 
ye ars with our da t aset. However, the photometric results 
of iGrankin et al.l ()2007l ). with observations secured over 
more than 20 years, find most Classical T Tauri stars to 
be fairly stable on long timescales. 

4. THE ORIGIN OF THE SCATTER IN ACCRETION RATE 
VS. MASS 

Our new constraints on accretion variability in a large 
sample allows us to constrain the origin of the distribu- 
tion of data points in the accretion rate vs. mass dia- 
gram. In Fig. m we show the two accretion indicators 
used in this study as a function of spectral type. For the 
Ha 10% width, we also include the non-accreting objects. 
In both cases, a clear trend is seen: the spectroscopic in- 
dicator drops towards later spectral types, which reflects 
the M vs. M correlation reported in previous studies (see 
As can be seen in the Ha 10% widths, the opposite 
trend is seen for the non-accretors, which is likely due to 
increasing levels of magnetic activity towards later spec- 
tral types. This leads to a 'U-shaped' distribution of 
data points: accretors and non-accretors are separated 
by a large gap at F~K spectral types, but are harder to 
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Fig. 3. — Change in Ca II-A8662 flux (upper panel) and 
H« 10% width (lower panel) vs. time. Data from Cha I are repre- 
sented by triangles, and those from Tau-Aur are drawn as squares. 
The maximum extent of the variability is reached after a few days 
and does not further increase towards longer timescales. 



distinguish in M-type objects. 

As can be seen in both panels, the total amount of scat- 
ter in the accretion indicator cannot be explained by in- 
trinsic accretion rate variability, which is shown as error 
bars. This is particularly obvious from the Ca II-A8662 
fluxes, which we identified as the more robust accretion 
indicator in Sj3l Converted to accretion rates, the ob- 
jects cover 1-2 orders of magnitude at any given spectral 
type, while the variability accounts only for a small frac- 
tion of this scatter (~ 0.35 dex, see SjS]). In other words, 
on timescales of days to months the objects do not move 
significantly in the diagram. 

Although we cannot exclude the presence of signifi- 
cant variability on longer timescales, the scatter in the 
M vs. M plot is most likely dominated by object-to- 
object varia t ions. This confirms previous claims by 
Natta et"an (120041) and is co ntrary to suggestions by 
Scholz &: Javawardhanal (120061) based on a much smaller 
sample. Thus, stellar mass or parameters that strictly 
scale with stellar mass are clearly not the only factors 
that affect the accretion rate. 

One possible explanation for the spread of accretion 
rates is the different evolutionary stages of the objects 
considered here. Accretion rates are expected to drop 
with age, either following a power law, if the timescale 
is determined by t he slow viscous evolution of the disk 
(jClarke fc Pringlel[2006f ) . or in a more rapid process, con- 
sistent with the rapid inner disk clearing inferred from 
the scarcity of 'transition' objects with optically thin 
inner disks. The presence of the clear gap between 
the accreting and non-accreting populations in Fig. 31 
particularly well-defined for the higher mass objects is 
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Fig. 4.— Ca II-A8662 flux (upper panel) and Ho 10% width 
(lower panel) vs. spectral types. Error bars indicate the full range 
of the variations. In the lower panel, we also show the non-accretors 
(hollow symbols) with Ha 10% widths < 200 km s"'^ (after sub- 
tracting rotational broadening). Objects without Ha emission are 
shown with a 10% width of km s~^. The r ight axis shows t he ac - 
cretion rates inferred from the relations fromlMohant v et al.l II2005I ) 
and Natta et al. (2004). For the M and K type stars, one sees a 
clear rise in Ha 10% width towards earlier spectral type, but at KO 
and earlier, it appears to plateau, with no objects showing widths 
larger than ~ 600 km s~^. 



thus highly interesting. The scarcity of objects in the 
transition from accretors to non-accretors argues for a 
rapid evolution between states, as suggested in disk evo- 
lution models inc l uding mass loss from the disk (e.g. 
IClarke et al.] 120011 : lArmitage et al.ll2003f ). This may be 
related to the similarly rapid timescale of transition be- 
tween objects with IR excess and those without. 

Apart from the evolutionary stage, a number of other 
factors have been suggested to influence the accretion 
rates. If disk ionization is a major factor, our findings 
imply that the main source of ionization is probably ex- 
ternal (e.g. cosmic rays), as the stellar ionising radiation 
is expected to scale with mass and thus would not allow 
for a wide range of accretion rates at constant mass. In 
general, our results favour scenarios where initial or en- 
vironmental conditions are more important than stellar 
parameters. For instance, a model where the dispersion 
in initial disk parameters, e.g. mass and radius, deter- 
mines the distribution of data points in the M vs. M plot , 
as recently described bv [Alexander fc Armitagd (|2006f l. 
would be consistent with our results. 

To constrain the infiuence of stellar and environmental 
parameters on the appearance of the M vs. M plot, it 
would be necessary to isolate them from evolutionary ef- 
fects. It would be of particular interest to investigate 
accretion rates versus other parameters along lines of 
constant mass in the M vs. M diagram. Unfortunately, 
for these tasks the currently available samples may not 
be sufficiently large. In principle, however, such analy- 
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ses have the potential to provide important information 
on the nature of the accretion process in yoimg stars. 
Finally, the causal relationships between our accretion 
indicators and accretion rates are not fully understood, 
and we use them merely as empirical tools, based on pub- 
lished correlations. Nevertheless, the fact that indicators 
such as Ca II-A8662 line flux correlate remarkably well 
with the current empirical framework of measuring ac- 
cretion rates gives some confidence that these indicators 
relate to physical reality. 
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